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ABSTRACT 

We examine the implications of Modified Newtonian Dynamics (MOND) on the large 
scale structure in a Friedmann-Robertson- Walker universe. We employ a "Jeans swin- 
dle" to write a MOND-type relationship between the fluctuations in the density and 
the gravitational force, g. In linear Newtonian theory, |g| decreases with time and 
eventually becomes < g , the threshold below which MOND is dominant. If the 
Newtonian initial density held has a power-law power-spectrum of index n < — 1, 
then MOND domination proceeds from small to large scale. At early times MOND 
tends to drive the density power-spectrum towards k" 1 , independent of its shape in 
the Newtonian regime. We use N-body simulations to solve the MOND equations of 
motion starting from initial conditions with a CDM power-spectrum. MOND with 
the standard value g = 10 _8 cm s -2 , yields a high clustering amplitude that can 
match the observed galaxy distribution only with strong (anti-) biasing. A value of 
g « 10~ 9 cms -2 , however, gives results similar to Newtonian dynamics and can be 
consistent with the observed large scale structure. 
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1 INTRODUCTION 

Perhaps the best evidence for the existence of dark matter 
is the flattening of rotation curves of many spiral galax- 
ies. Newton's g N ~ l/(distance) 2 law for the force of grav- 
ity predicts that the observed distribution of baryonic mat- 
ter should show a fall off in the rotation curves, and hence 
the need for dark matter to explain the flattening. A uni- 
verse predominantly filled with cold dark matter (CDM) is 
consistent with observations of the cosmic microwave back- 
ground, and the galaxy population. Nevertheless, the dark 
matter scenario does not readily explain the shape of rota- 
tion curves of some galaxies in the inner regions (e.g., Moore 
1994, Flores & Primack 1994, de Blok et al. 2001), the dark 
matter particle has not been discovered, and we lack a di- 
rect experimental verification of Newton's law of gravity at 
low accelerations. So it seems prudent to examine alterna- 
tive scenarios explaining the flattening of rotation curves. 
Modified Newtonian Dynamics (MOND) acting at low accel- 
erations (e.g., Milgrom 1983, Bekenstein & Milgrom 1984) 
can explain the flattening of rotation curves without the 
need for dark matter. Two main versions of MOND exist. 
The first replaces Newton's law of gravity at g N < g g ^ by 
9 ~ \/9o9n where g « 1.2 x 10~ 8 cm s -2 is found to give the 



In MOND literature a is often used instead of g . 



best results in fitting the rotation curves. The second mod- 
ifies the law of inertia at low accelerations. Here we adopt 
the first version where the equations of motion can be de- 
rived from a Lagrangian, admitting conservation of energy 
and momentum (Bekenstein & Milgrom 1984). Because of 
the success of MOND at fitting the rotation curves (e.g., 
Milgrom & Braun 1988, Begeman et al. 1991, Sanders 1996, 
de Blok & McGaugh 1998, Sanders & Verheijen 1998), it 
is worthwhile to push it further and confront it with obser- 
vations of the large scale structure in the Universe (Felten 
1984, McGaugh 1999, Sanders 2001). A problem here arises 
however. MOND, when applied to a uniform background, 
predicts the collapse of any finite region in the universe 
regardless of the mean density in that region (e.g., Felten 
1984, Sanders 1998). To solve this problem Sanders (2001) 
proposed a two-field Lagrangian based theory of MOND in 
which the Friedmann-Robertson- Walker (FRW) background 
cosmology remains intact in the absence of fluctuations. He 
argued that this theory leads to large scale structure resem- 
bling Newtonian dynamics with CDM-like initial conditions. 
Here we take a simpler approach in which we employ the 
Jeans swindle to write a MOND type relationship between 
the fluctuations in the density and the gravitational force 
field. Although our procedure can be derived from a La- 
grangian and is equivalent to Sanders' two-field theory in 
the limit of small coupling, we do not see any physical jus- 
tification for it. Nevertheless, the equations of motion have 
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the desired properties and we proceed to explore their con- 
sequences on the evolution of large scale structure. 



2 MODIFICATION OF THE COSMOLOGICAL 
NEWTONIAN EQUATIONS OF MOTION 

The background FRW cosmology is described by the scale 
factor a(t) normalized to unity at the present, the Hubble 
function H(t) = a/a, and the total mean background den- 
sity ptot — p + A/(8irG), where p is the mean matter den- 
sity and A is the cosmological constant. Further, we define 
£1 — p/pc, where p c = 3H 2 /(8ivG) is the critical density. 
These cosmological quantities are related by Einstein equa- 
tions of general relativity. Let r and x = r/a denote, respec- 
tively, physical and comoving coordinates. The fluctuations 
over the uniform background in the matter distributions are 
described by the the comoving peculiar velocity u = dx/di 
of a patch of matter, the density contrast 5(x) = p(x)/p— 1, 
where p(x) is the local density, and the fluctuations in the 
gravitational force field, g. Neglecting thermal effects, the 
Newtonian equations of motion governing the evolution of 
the fluctuations are: the continuity equation 

^+V„.(l+*)u = 0, (1) 



the Euler equation of motion, 



du 

dt 



+ 2Hu = g/a , 



and the Poisson equation, 



1 



V x • g = -4ttG/9(5 



-SIH J 5 . 



(2) 



(3) 



In the linear regime these equations imply that g| decreases 
with time (see 3.1). Once |g| drops below g , MOND takes 
over the subsequent evolution of the fluctuations. The con- 
tinuity equation (jl|) reflects mass conservations and holds 
also in MOND. The MONdification (terminology by M. Mil- 
grom) must be done by either changing the law of inertia, 
i.e., the Euler equation (^), or the relation between the den- 
sity and the force field, i.e., the Poisson equation Here 
we choose to maintain the Euler equation and replace the 
Poisson equation with 



-QH 2 5 



(4) 



In writing this equation we have assumed that MOND only 
affects the fluctuations and leaves the background cosmology 
intact. This "Jeans Swindle" (Binney & Tremaine 1987) is 
hard to justify on physical grounds, even though it leads to 
equations of motion which can be derived from a Lagrangian, 
admitting conservation of energy and momentum (Sanders 
2001, Bekenstein & Milgrom 1984). Although we can take 
g to be a function of time here we assume it is a constant 
and write it as 



9o = f cH o = fRaH a , 



(5) 



where / is a numerical fudge factor, c is the speed of light, H 
is the current value of the Hubble function, and Rh = c/H 
is the Hubble radius. Fitting the rotation curves by MOND 
requires / « 1/6 for H = 70km s _1 Mpc -1 (e.g. Sanders 
1996), but we will see that / must be smaller for MOND 



to be in reasonable agreement with the observed large scale 
structure. Equation mft determines g only up to a divergence 
free vector field h. So in a sense the theory is incomplete. 
However, Bekenstein & Milgrom (1984) have shown that h 
decays rapidly with increasing scale. So on large scales we 
neglect h and assume that jgjg is irrotational so that it is 
fully determined by (^J). For conciseness we gather all the 
time dependent function in the Euler equation. So let us 
define an irrotational vector field g N satisfying 



V • g N = -S 
Then, by @ 
'3 



we have 



aQH g 



1/2 



,2— ~ J V |g N |i/ 2 ' 
So the Euler equation becomes 



du 



\ 2a 



1/2 



\S, 



11/2 



The Newtonian equivalent of this equation is 



du 

m 



+ 2Hu 



(6) 



(7) 



(8) 



(9) 



3 THE LIMIT OF SMALL DENSITY 

FLUCTUATIONS AND EARLY TIME MOND 

We treat now the MOND equations in the limit of small 
density fluctuations. We focus on the evolution at early 
times when Q « 1, restricting the treatment to a matter 
dominated universe. Assuming that the fluctuations are ini- 
tially Newtonian, we discuss the epoch at which they become 
MOND dominated, and the shape of their power-spectrum 
in the MOND regime. 

The modified Poisson equation (Q) relating g and 8 
is non-linear and cannot be linearized. However, we can 
write the equations at early times when S and|V • u/H\ 
are still C 1. In this limit we neglect the term V • (<5u) 
in the continuity equation (|l]) and replace the full derivative 
du/ dt = du/ dt + u • Vu with du/ dt in the Euler equation 
(^). The continuity equation then implies that 8 — —V ■ u 
which in combination with equation (p|) gives 



dt ' 

Substituting this result in the Euler equation 
glecting u • Vu in du/ dt yield 



(10) 



and ne- 



d 2 s 



at'' at 



1/2 



I Si- 



ll/2 



(11) 



This system of 3 first order differential equations determines 
g N locally at any point in space. The full solution depends 
on the values of g N and its time derivative at some ini- 
tial time, ti. We could not find analytic solutions to this 
equation, but in an Q = 1 universe the asymptotic behavior 
t^$>U can easily be obtained by substituting g = b t a . This 
gives 

isj = ^9 tytt) , (12) 



where t is the Hubble time. Using ([?]), this implies 



|g| = Jq9 



(13) 
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redshift at which fluctuations on a comoving scale x become 
MOND dominated. Assume that the transition epoch be- 
tween Newtonian dynamics and MOND occurs early enough 
so that Q w 1, an assumption that will be justified below. 
The initial density field in the Newtonian regime is gaussian 
and characterized by the power-spectrum, < \8k\ 2 >, where 
5k is the Fourier transform of the density as a function of 
the wavenumber, k. If the power-spectrum is a power-law of 
index —3 < n < 1, then the rms value of the density field 
smoothed on a comoving scale x at redshift z is approxi- 
mately 



<t(x, z) — a(z) I — 



(14) 



where the redshift dependence through a(z) — 1/(1 + z) 
follows from linear Newtonian theory if fl = 1, and the scale 
x„ sets the amplitude of the fluctuations. According to either 
(H) or g N ~ GSM/r 2 , the rms of g N on a scale x is 



-H axa 



(15) 



Figure 1. The absolute value of g/g in the limit of small density 
fluctuations as a function of a in an Q = 1 matter dominated 
universe. The solid curve is the behavior according to Newtonian 
dynamics. The dotted curve is obtained when MOND is activated 
at |g| < So- 



So |g| approaches a constant that is smaller than g , 
i.e., fluctuations entering MOND never become Newtonian 
again. Integration of the MOND equations for the spherical 
collapse model (not shown here), and the N-body simula- 



tions presented in section 4 show that|g| < g also in the 
moderately non-linear regime. In Fig. hi we plot |g|/<? D as a 
function of a, in the limit of small density fluctuations in an 
f2 = 1 universe. The solid line follows the growing mode of 
linear Newtonian theory. The dotted line shows |g| in the 
MOND regime and is obtained by numerical integration of 
( pd| ) for |g| < g . This numerical solution confirms our an- 
alytic result that |g| approaches (3/10)/? . This asymptotic 
behavior has interesting implications on the topology of the 
density field. It implies that one-dimensional pancake-like 
perturbations the density contrast 8 = —V ■ g N vanishes 
in all space except at the surfaces of g = 0, where it can 
be described by a Dirac-delta function. In a generic three- 
dimensional perturbation, surfaces of g = are still singular 
but the density contrast in the rest of space is non-vanishing 



and grows with time like 8 ~ t 



4/3 



as compared to <5 ■ 



for the growing mode in linear Newtonian theory (e.g. Pee- 
bles 1980). We will use this result in 3.2 to estimate the 
shape of the power-spectrum in early time MOND. 



3.1 Epoch of MOND domination 

The fluctuations in the Newtonian gravitational force field 
over a spherical region of physical radius r = ax and density 
contrast 8 are g N ~ G8M /(ax) 2 , where M ~ p(ax) 3 . Since 
p ~ 1/a 3 , g N decreases with time if 8 ~ t a with a < 2. 
In linear Newtonian theory, a < 1 and so at early enough 
times g N > g and the fluctuations would be in the Newto- 
nian regime. At later times g N drops below g and the fluc- 
tuations enter the MOND regime. Here we determine the 



The redshift, z mond , at which g N = g = fRnH 2 is then 
:.. ■ = - J I — I 2 -1 ■ (16) 



2 to x 

o Xq \ X, 



So if n < —1 then fluctuations on small scales enter the 
MOND regime before large scale fluctuations. For n > — 1 
the interpretation that large scales become MOND domi- 
nated earlier is somewhat vague. It is hard to visualize a 
large scale fluctuation in the MOND regime while the su- 
perimposed small scale fluctuations are still Newtonian. A 
more appropriate interpretation is that for n > — 1 the vol- 
ume of space containing fluctuations in the MOND regime 
is mainly made up of large regions separated by small scale 
fluctuations that are still Newtonian. This interpretation is 
sustained by the following argument. Since in the Newto- 
nian regime g is gaussian, the probability that a spherical 
region of radius x becomes MOND dominated at time t is, 



P oc 



exp i 



-y 2 /2) y 2 dy 



According to a s 



-(n + l)/2 



(17) 



so if n > —1 then P 



increases with scale, meaning that the fraction of space al- 
ready in MOND is made up of large regions. If n < — 1 then 
P decreases with scale and MOND penetrates small regions 
before it dominates the whole space. 

Fluctuations enter the MOND regime at a relatively 
early time. Taking i„~ lOh^Mpc and R H ~ 3000h _1 Mpc 
then, according to (Jlq) , the transition into MOND occurs 
at z > /200, on scales ^ 10h _1 Mpc if n > —1, and 
10h - Mpc if n < —1. In a CDM-like power-spectrum the 
effective power-index on scales ^ 20h _1 Mpc is n ^ — 1 and 
changes to —3 < n ;$ — 1 on smaller scales. For this type of 
power-spectrum the transition into MOND of all scales oc- 
curs at an early time when the density fluctuations are still 
small and £1 ~ 1. This considerably simplifies our task of es- 
timating the density power-spectrum in the MOND regime 
in the next subsection. 



© 0000 RAS, MNRAS 000, 000-000 



4 Nusser 



3.2 power-spectrum in the MOND regime at 
early times 

MOND gives accelerations oc ^/g Q g N > g N , leading to 
faster growth rates than Newtonian gravity. So we expect 
a break in the shape of the power-spectrum at the scale just 
entering MOND. Assuming a power-law power-spectrum, 
< | Sk\ 2 >~ k n , in the Newtonian regime we estimate the 
power-index of the density fluctuations after they become 
MOND dominated. Consider the rms value, a(x, a), of the 
density field over any scale x. Suppose that at redshift 
2i > 1 fluctuations on all scales were Newtonian so that 
a — a(x,Zi), where a is given by (|l4|). Let x maiid (z) be the 
scale just entering the MOND regime at redshift z < Zi. For 
simplicity we consider the case n < — 1, so at x > x ond 
the evolution is Newtonian and a = a. A fluctuation on 
a scale x < x mond entered MOND at redshift z x > z. As 
we have seen in the previous subsection fluctuations on a 
wide range of scales become MOND dominated at an early 
enough time so that fi~l.Inanf2~l universe, the growth 
factors in the Newtonian and MOND regimes are 8 ~ a and 
5 ~ a 2 , respectively. So, from z; until z x the fluctuation 
grew by a factor of {a x /a%), while from z x until z it grew by 
(a/a x ) 2 , where a% = a(zi), a x — a(z x ), and a = a(z). So for 
x < x ,, 



a(x,a)=(— ) (—\o(x, ai ) 
\a x J \a,i J 

According to ( |l6|) 
a x \ x J 



and 



rrl x, a) = a 



(18) 



(19) 



(20) 



so the power- index in the MOND regime is n mond = — 1. 
A similar argument for n > —1 yields the same result. 
Therefore the power-index in the MOND regime is inde- 
pendent of the value of n. A general power-spectrum can 
be approximated as a power-law over any limited range of 
scales where our result applies. Therefore, MOND drives the 
power-spectrum towards independent of its shape in the 
Newtonian regime. 

The treatment here ignored mode-coupling and as- 
sumed that all fluctuations became MOND dominated when 
fl was still very close to unity. Although the simulations pre- 
sented in the next section do not show strong mode-coupling 
for low values of g , the shape of the power-spectrum derived 
here should serve as a general indication only. 



4 RESULTS FROM N-BODY SIMULATIONS 

In this section we study the general evolution of fluctua- 
tions under MOND by means of N-body simulations. This is 
done by adapting a particle-mesh (PM) N-body code orig- 
inally written by E. Bertschinger to solve the Newtonian 
equations (Bertschinger & Gelb 1991). At each time step, 
the code readily provides the irrotational component of g N 
from the density field as estimated on a grid from the parti- 
cle distribution. We adapt the code by adding a routine to 



Table 1. Details of the six simulations. Each simulation contained 
64 3 particles in a cubic box of 128h _1 Mpc on the side. The rms 
value, <rg, of density fluctuations in spheres of radius 8h Mpc is 
measured from the output of each simulation at the final time. 

simulation gravity g /(1.2 X 10 — 8 cm/s 2 ) f2o erg 



I 


MOND 


1 


1 


3.9 


II 


MOND 


1 


0.03 


4.83 


III 


MOND 


1 

(10- 6 5 initial ) 


1 


1.16 


IV 


MOND 


1/12 


1 


1.08 


V 


MOND 


1/12 


0.03 


1.01 


VI 


Newton 





1 


0.86 



Figure 2. The particle distribution in a slice of thickness 
2.56h~ 1 Mpc through the simulations at the final time. The pan- 
els in a clockwise order correspond to simulations I — VI listed in 
the table 1. See flg2.gif 



compute the gravitational field g from g N using the relation 
([])• We also incorporate the time dependent functions in the 
MOND equations in a leapfrog time integration scheme to 
move the particles to the next time step. 

Five MOND simulations with different values of g a 
and fio = = 1) were run. For comparison, a simula- 
tion with Newtonian gravity was also run. The simulations 
started with gaussian initial conditions generated from the 
CDM power-spectrum with H = 50h _1 Mpc and flo = 1. 
Each simulation contained 64 3 particles in a cubic box of 
128h _1 Mpc on the side. The CDM power spectrum on the 
large scales probed by the box is close to k~ l as motivated 
by our result for the density power-spectrum in the MOND 
regime at early times. The parameters of each simulation 
are written in Table 1. Except simulation III, all simula- 
tions started with the same initial conditions which were 
normalized so that the density field in Newtonian simula- 
tion at the final time had an rms value of as = 0.86 in 
spheres of radius 8h~ Mpc. In simulation III, the amplitude 
of the initial conditions was reduced by a factor of 10 6 . Sim- 
ulations with flo = 1 and 0.03 started at a = 0.049 and 
0.01, respectively. At these two time the linear Newtonian 
density growth factor, D, was the same for the two values 
of ilo ■ For the high g , this choice of the starting time of the 
simulations corresponds roughly to the epoch at which fluc- 
tuations on all scales become MOND dominated according 
to our estimate in section 3.2. The particle distribution at 
this time was obtained using the Zel'dovich approximation. 
We are interested in comparing the net dynamical effect of 
decreasing g and so we start MOND simulations with low 
and high g at the same value of D. 

Fig. ^| show the particle distribution in a slice in the 
simulations at a = 1. This figure demonstrates that MOND 
with the standard value of g ~ 1.2 x 10 _8 cms -2 (bottom 
and middle panels to the left) yields high clustering at the 
present time, for both fio = 1 and 0.03. Lowering the ampli- 
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Figure 3. The variance A 2 (fc) = k z < \5 k \ 2 > /{2ir 2 ) of den- 
sity fields per ln(fc) from the simulations at the final time. The 
solid, long-dashed, short-dashed, dotted, dot-short-dashed, and 
dot-long-dashed and lines correspond respectively to simulations 
I-VI. 



Figure 4. The peculiar velocity vs. the "Newtonian" gravity field 
g N obtained from the density contrast using (^) . The panels show 
fields from simulations If (Oo = 0.03, high g ), fV (firj = f, low 
g ), V (f2o = 0.03, low g ), and VI (Newtonian), as indicated 
in the figure. The fields were smoothed with a Gaussian window 
of radius 8h — iJVlpc. The diagonal line in each panel is plotted to 
guide the eye. See fig4.gif 

tude of the initial fluctuations by a factor of 10 reduces the 
clustering at the present time but severely washes out small 
scale fluctuations (top panel to the left). For both values 
of Slo, MOND with a smaller value for g produces particle 
distribution (top and middle panels to the right) similar to 
Newtonian dynamics (bottom panel to the right). The values 
of as measured from simulations at the final time are listed 
in Table 1. Fig. | shows A 2 (fc) = Vk 3 < \S k \ 2 > /2n 2 where 
V is the volume of the simulation box (Peacock 1999). This 
quantity gives the contribution to the variance per ln(fc). 
The MOND curves with low g (IV and V) follow the shape 
of the Newtonian curve (VI). 

A common method for estimating the cosmological den- 
sity parameter, f2, is the comparison between the measured 
peculiar velocity field and the "Newtonian" gravity field, 
g N , inferred observed distribution of galaxies using (^) (see 
Strauss & Willick 1995 and references therein). Putting aside 
the issue of redshift distortions, linear Newtonian theory 
gives /(fi)u = Hg N , where /(fi) « fi 6 (e.g., Peebles 1980). 
This comparison has been done using several data sets with 
the conclusion that the observed velocity and gravity fields 
are tightly correlated. To test whether MOND can account 
for the correlation seen in the observations, we have per- 
formed a comparison between the velocity and "Newtonian" 
gravity field, g N , computed from the simulations at the final 
time. The results are shown in Fig. 0. To avoid strong non- 



linearities, the fields were smoothed with a gaussian window 
of radius 8h _1 Mpc. The rms values of the one dimensional 
u and Hg N shown in the figure are (a u ,a g ) = (1176, 336), 
(450,302), (223,336), and (209,226) in kms" 1 , respectively, 
for the simulations II, IV, V, and VI. Simulation II is highly 
non- linear even on scales larger than 8h _1 Mpc and the rel- 
ative scatter between the fields is large. In the simulations 
with small g , the correlation between the fields is tight for 
both fio = 1 (simulation IV) and flo = 0.03 (V). The scatter 
between the fields is larger than in the Newtonian simulation 
(VI), but it is negligible compared to observational uncer- 
tainties. The slope of the regression of u x on Hg N x is higher 
in simulation IV than in the Newtonian simulation, VI. 
These two simulations have tto = 1 but as is the case with 
rotational speeds of galaxies, MOND also accounts for larger 
velocities on large scales. These results indicate that MOND 
can be consistent with this type of analysis. In MOND, the 
slope of the regression of u on Hg N is also a decreasing 
function of fi, but is not proportional to /(fi). For example 
for the MOND simulation V with low g and Qq = 0.03, 
the slope is 0.55, instead of 0.25 in linear Newtonian the- 
ory. Note that the relation (Q) between g and g N implies 
that |g| go in all MOND simulations, sustaining our con- 
jecture that large scale fluctuations never leave the MOND 
regime. We have also computed the probability distribution 
function (PDF) of density field smoothed on 8h _1 Mpc. The 
shape of the PDF in the simulations with small g is very 
close to that in the Newtonian simulation. 

The simulations shown so far were run assuming all par- 
ticles enter the MOND regime at the same time. This cap- 
tures the net dynamical effects of varying g . However, a 
realistic MOND simulation should start at an early enough 
time when most particles are still Newtonian. A particle 
moves accroding to Newtonian dynamics, until the its pecu- 
liar gravitational (in absolute value), g, drops below g . The 
subsequent motion of the particle is governed by MOND as 
long as g < go- We have run a suite of simulations with this 
recipe for modeling the transition into MOND. Of course the 
transition into MOND could be made a little smoother but 
in the lack of an exact model for how the transition occurs 
we choose to proceed assuming sudden transition. The simu- 
lations we have run correspond to the CDM power-spectrum 
used in the simulations shown in Fig. ^| and two scale-free 
power-spectra with power-indices, n = —1 and n = 1. All 
simulations were run for Q — 0.03 with the standard and 
low value of go- The initial conditions were normalized sim- 
ilarly to the Newtonian simulation in table 1. Slices of the 
particle distribution in these simulations are shown in Figs. 
|H] and ^| for a = 1 and a = 0.1, respectively. The particle 
distribution of the CDM simulation in Fig. M (bottom pan- 
els) is a little more evolved than in Fig. ^. Despite this, the 
particle distributions in the two simulations are very similar. 
This indicates that MOND becomes dominated at a similar 
epoch for all particles and that our estimate of this epoch 
used in the simulations shown in Fig. B is reliable. 



5 SUMMARY AND DISCUSSION 

Several arguments against MOND have been presented in 
the literature ( e.g., Felten 1984, van den Bosch et al. 2000, 
Mortlock & Turner 2001, Aguirre et al. 2001, Scott et al. 
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Figure 5. The particle distribution in a slice of thickness 
2.56h _1 Mpc through the simulations at the final time, a = 1. 
These simulations have = 0.03 and started early enough when 
most particle were Newtonian. A particle became MOND domi- 
nated once its peculiar acceleration dropped below go . The panels 
on the left are for simulations run with standard go while the pan- 
els on the right are with low go (1/12 standard). Bottom, middle, 
and top rows correspond to power-spectra of CDM, n = —1, and 
n = +1, respectively. See flg5.gif 



Figure 6. The same as the previous figure 
at a = 0.1. See fig6.gif 



, but for simulations 



2001). The reader may decide for her/himself which of these 
arguments are compelling. Here we focus on examining the 
general properties of large scale structure under MOND 
and put aside the arguments presented in the literature 
against MOND. The value of g necessary to fit the rota- 
tion curves of galaxies produces far too much clustering and 
can be consistent with observations only with strong anti- 
bias. Decreasing g has two effects: 1) it pushes the epoch of 
MOND domination to later time, and 2) it reduces the typ- 
ical accelerations produced by MOND forces. The simula- 
tions presented in this paper show that MOND with smaller 
g leads to large scale structure similar to Newtonian dy- 
namics. This is encouraging for MOND since the standard 
paradigm based on Newtonian gravity and dark matter is 
successful at explaining many observations of the large scale 
structure and the galaxy population (e.g., Kauffman, Nusser 
& Steinmetz 1997, Benson et al. 2000, Diaferio et al. 2001) 
. Still one has to adopt a generalization of MOND, allowing 
for, or at least mimicking, a scale dependent g . Sanders' 
(2001) two-field theory for example can serve as a general 
framework for that purpose. In this theory the rapid growth 
of fluctuations in MOND is tamed by the coupling between 
the two fields. The amplitude of the initial conditions used 
in the simulations roughly matches the normalization im- 
plied by measurements of the cosmic microwave background 
anisotropies (Hu W., Sugiyama & Silk 1997). We have seen 
that by reducing the amplitude of the initial conditions by a 
large factor MOND with the standard value of g produces 
the correct as, but causes a severe washing out of structure 
on scale of 10s of h _1 Mpc. Therefore it seems necessary to 
decrease g to match the large scale structure. 

Our goal here was to present general criteria for the con- 
sistency of MOND with large scale structure. Further study 
of MOND within a realsitic recipe for galaxy formation to 
clarify the role of (anti-) biasing remains to be done. 
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